Raman spectroscopy of the interlayer shear mode in few-layer M0S2 flakes 
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Single- and few-layer M0S2 has recently gained attention as an interesting new material system 
for opto-electronics. Here, we report on scanning Raman measurements on few-layer M0S2 flakes 
prepared by exfoliation. We observe a Raman mode corresponding to a rigid shearing oscillation 
of adjacent layers. This mode appears at very low Raman shifts between 20 and 30 cm -1 . Its 
position strongly depends on the number of layers, which we independently determine using AFM 
measurements and investigation of the other characteristic Raman modes. Raman spectroscopy of 
the shear mode therefore is a useful tool to determine the number of layers for few-layer M0S2 flakes. 



The tremendous growth of experimental research on 
graphene in the past few years stems, in part, from 
the simple exfoliation technique that allows for prepa- 
ration of single- and few-layer flakes from bulk crys- 
tals. This technique is applicable to many layered crystal 
structures in which the binding energy between adjacent 
planes is much lower than the binding energy within a 
planed. Among these layered structures, the dichalco- 
genide M0S2, which is used commercially, e.g., as a high- 
temperature dry lubricant, has attracted a lot of interest. 
It was recently shown to undergo a transition from indi- 
rect to direct-gap semiconductor when its thickness is 
reduced to a single layes 2 ^, leading to pronounced pho- 
toluminescence. This drastic change of the band struc- 
ture has been investigated theoretically by a number of 
groups^—, and further calculations suggest the possibil- 
ity of band structure engineering using strain 7 -. Low- 
temperature photoluminescence measurements revealed 
the presence of impurity-bound excitons in single-layer 
M0S2 flakes 8 -, which can be suppressed in oxide-covered 
M0S2 9 . The photocarricr lifetime in single-layer M0S2 
is sufficiently short 8 - to make the material interesting for 
fast photodetectors, and a MoS2-based phototransistor 
was reported recently^. Room-temperature transistor 
operation with very large on/off ratio has been demon- 
strated for single M0S2 layers^. As in graphene, where 
research was initially focused on single layers and later 
expanded to study also bilayers and trilayers due to their 
different band structure, few-layer M0S2 flakes may be 
interesting for their transport properties: very recently, 
ambipolar transistor operation was shown in few-layer 
M0S2 12 . Both, chemical exfoliation-^ and vapor phase 
growth techniques-^ have been demonstrated for M0S2, 
indicating the possibility of fabricating large-area thin 
films neccessary for potential applications. Similar to 
graphene, where Raman scattering has been used to de- 
termine, e.g., the layer thickness^ or the doping type and 
concentration 16 , Raman spectroscopy is a highly useful 
tool to map M0S2 flakes and to identify single layer o 17 ' 18 . 
This is facilitated by the fact that two characteristic Ra- 
man modes, A± g (out-of-plane optical vibration of the 
sulfur atoms) and (in-plane optical vibration of Mo 
and S atoms), show opposite frequency dependence on 



the number of layers: the Ai g mode increases its fre- 
quency with the flake thickness, while the E^ anoma- 
lously softens due to increased dielectric screening^, so 
that the difference of the mode frequencies is character- 
istic for a certain number of layers. While many of the 
Raman modes of graphene and M0S2, or other dichalco- 
genides, are very different due to the different crystal 
structures, two modes are rather generic for layered crys- 
tal structures: 

(1) a shear mode, in which adjacent layers rigidly oscil- 
late relative to each other, with the oscillation amplitude 
lying in the layer plane. This mode was observed for, 
both, bulk graphite 20 , and bulk MoS^- More recently 
it was detected in few-layer graphene 2 ^. 

(2) a compression mode, in which adjacent layers oscil- 
late rigidly, with the oscillation amplitude perpendicular 
to the layer plane. This mode was recently observed in 
few-layer graphene in a double resonant process in com- 
bination with an LO phonon 2 ^. 

Here, we report on scanning Raman scattering mea- 
surements on M0S2 flakes. The M0S2 flakes were pre- 
pared using the transparent tape liftoff method well- 
established for graphene, from natural M0S2. A p-doped 
silicon wafer (specific resistivity p — 0.005 Ohm cm, cor- 
responding to degenerate doping) with 300 nm SiC>2 layer 
and lithographically defined metal markers was used as 
a substrate. After initial characterization with an opti- 
cal microscope, the samples were analyzed by Raman 
spectroscopy at room temperature. For this, we uti- 
lized a microscope setup, in which a 532 nm cw laser 
was coupled into a lOOx microscope objective, which also 
collected the scattered light in backscattering geometry. 
The scattered light was recorded using a triple grat- 
ing spectrometer equipped with a liquid-nitrogen-cooled 
charge-coupled device (CCD) sensor. The sample was 
mounted on a piezo-stepper table and scanned under the 
microscope. The spatial resolution of this setup is about 
500 nm. For some of the Raman measurements, a cross- 
polarized backscattering geometry was used. Here, a po- 
larizer was placed in front of the spectrometer so that 
only the scattered light with polarization perpendicular 
to the laser polarization was coupled into the spectro- 
meter. In this scattering geometry, the spectrally broad 
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FIG. 1: (a) Raman spectra of a three-layer M0S2 flake. The 
top spectrum is recorded without a polarizer in front of the 
spectrometer, the bottom spectrum in cross-polarized geom- 
etry, (b) Normalized Raman spectra of the shear mode mea- 
sured on M0S2 flakes of different thickness, (c) Spectral posi- 
tion of the shear mode as a function of (number of layers) - . 
Multiple datapoints indicate measurements on different flakes. 
The dotted line indicates the mode position as predicted by 
the model developed in ref. (2^ 
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FIG. 2: (a) Optical micrograph of a few-layer M0S2 flake, 
(b) AFM image of the same flake. The two lines indicate the 
areas where the height profiles (shown in (c)) were taken, the 
dots indicate the points where the Raman spectra (shown in 
(d)) were recorded, (c) AFM height profiles for steps from 
the substrate onto different areas of the flake, (d) Raman 
spectra measured at two different points on the flake. The 
dotted lines serve as guide to the eye. 



background, which stems from inelastic light scattering 
of free carriers in the heavily p-doped Si substrate^, is 
suppressed and allows us to obtain Raman spectra very 
close to the laser line. 

First, we discuss the typical Raman spectra observed 
in our samples. Figure Q] (a) shows two Raman spec- 
tra collected on a three-layer area of a flake. The top 
spectrum is recorded without a polarizer in front of the 
spectrometer, while the bottom spectrum is recorded in 
cross-polarized geometry, as described above. Due to the 
large and spectrally broad background caused by the sil- 
icon free-carrier scattering, the unpolarized spectrum is 
only recorded in a spectral range starting at 200 cm -1 
with respect to the laser line, while the cross-polarized 
spectrum is collected in a range starting at 10 cm -1 . In 
the unpolarized spectrum, we note the two characteristic 
Raman modes (386.1 cm -1 ) and A lg (409.7 cm" 1 ), 
with a frequency difference of 23.6 cm -1 as expected for 
three layers^. By contrast, in the cross-polarized spec- 
trum, the Ai s mode is almost completely suppressed, 
while the Yj\ g mode remains visible. Additionally, in 
the cross-polarized spectrum, we observe a low-frequency 
mode at 25.3 cm -1 , which we associate with the E| 9 in- 
terlayer shear mode. The main focus of this manuscript 
is the study of its frequency dependence on the number 
of M0S2 layers. 

Our results are summarized in Figure [If b): it shows 
normalized shear mode spectra obtained from M0S2 
flakes of different thickness. The number of layers was 
determined independently, by, both, atomic force micro- 



scopy (AFM) measurements (see below), and investiga- 
tion of the Ai g - frequency difference. We clearly ob- 
serve that the shear mode shifts from below 20 cm -1 for a 
bilayer to about 30 cm -1 for bulk-like flakes, a very large 
relative shift of more than 50 percent. In Figure [He), we 
plot the shear mode positions extracted for 16 different 
flake areas as a function of (number of layers) -1 . Multi- 
ple data points for a given thickness stem from measure- 
ments on different flakes. For comparison, we have calcu- 
lated the expected mode position as a function of (num- 
ber of layers) -1 as predicted in ref . [22] using a linear chain 
model with fixed interlayer coupling strength. Here, we 
have used our measured value for the shear mode fre- 
quency in bulk as the only parameter. We clearly see 
that our experimental values show a more pronounced 
shift than expected using the linear chain model. Two 
factors may contribute to this difference: (a) we have to 
assume that, since all experiments were carried out un- 
der ambient conditions, the topmost layer of the M0S2 
flake is partially covered with adsorbates, effectively in- 
creasing the mass per unit area, (b) our flakes are not 
free-standing, but deposited on a SiC>2 substrate, there- 
fore, the bottom layer will experience a weak van der 
Waals interaction with the substrate, which would act as 
an additional spring constant in the linear chain model. 

Figure [2] demonstrates the process of establishing the 
relation between shear mode position and number of lay- 
ers in more detail. Utilizing flakes with large areas of 
uniform height, determined by optical microscopy, such 
as the one shown in Fig. [Ha) , we determine the thick- 
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FIG. 3: (a) Optical micrograph of a few-layer M0S2 flake. 
The square marks the scan region for the false color plots in 
(b) and (c). (b) False color plot of the frequency difference 
between the Ai 9 and E2 9 modes. The outline of the single- 
layer area of the flake is marked by the solid line, (c) False 
color plot of the shear mode position. The scan region is 
identical to (b). 



ness of these areas by extracting height traces from AFM 
images. Raman spectra are then recorded at points on 
the flake for which the height is known. In Fig. &b) , the 
AFM image corresponding to Fig. HJa) is shown. The 
two colored lines indicate the position of the AFM line 
traces leading from the substrate to the flake shown in 
Fig. [H[c). From these traces, using the thickness of a sin- 
gle M0S2 layer of ^0.7 nm, we find that the black trace 



corresponds to a step height of 5 layers, while the red 
trace corresponds to 4 layers. Raman spectra collected 
on these areas of the flake (Fig. Hfd)) demonstrate the 
large blueshift of the shear mode (about 1.8 cm -1 ), which 
is clearly visible in the data and also easy to determine 
using an automated fitting routine. By comparison, the 
F,2 g mode only shows a weak redshift (about 0.5 cm -1 ), 
so that even though the spectrally integrated intensity 
of this mode is about 8 times higher than that of the 
shear mode for a 5-layer region, the shear mode allows 
for easier discrimination of the number of layers. 

The large relative shift of the interlayer shear mode 
can be used to efficiently map few-layer M0S2 flakes by 
using the shear mode position. One example is shown 
in Figure [3] The few-layer flake shown in Fig. [3ja) was 
scanned in the region marked by the square using the 
Raman microscope with a step size of 300 nm. The 
collected spectra were analyzed using an automated fit 
routine, which determines the peak positions of the Ai g , 
E^g and E| s modes as a function of position. From this 
data, false color plots were generated. In Fig. [^b) , we 
plot the frequency difference Ai g - E3, 9 as a function of 
position. From the frequency difference we can clearly 
identify regions of different thickness, including a small 
monolayer area, which shows a frequency difference of 
about 18 cm -1 and is marked with a solid outline. Na- 
turally, there is no shear mode signal in this region in 
Fig. EIc), where we plot the shear mode frequency as a 
function of position. Comparing the contrast of the two 
false color images, where we have used the same false 
color scale, we immediately notice that the shear mode 
image shows much more well-defined transitions between 
regions of different thickness and easily allows to iden- 
tify, e.g., a small trilayer region (green color) within the 
larger bilayer area (blue color), and a transition from 
6 layers (dark yellow color) to about 10 layers (orange 
color) , which are hard to discern from the frequency dif- 
ference image. 

In conclusion, we have used scanning Raman spec- 
troscopy to study the interlayer shear mode in few-layer 
M0S2 flakes. We observe a large frequency shift of the 
shear mode to higher energies with increasing number of 
layers. This large shift allows us to precisely map the 
layer thickness in few-layer M0S2 flakes. The authors 
gratefully acknowledge financial support by the DFG via 
SFB689, SPP 1285 and GRK 1570 as well as fruitful dis- 
cussion with L. Wirtz, A. Molina-Sanchez, J. Maultsch 
and N. Scheuschner. 
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